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Abstract 

The transition from arbitrary to chaotic fluctuation properties in quantum sys- 
tems is studied in a random matrix model. It is assumed that the Hamiltonian 
0^ can be written as the sum of an arbitrary and a chaos producing part. The 

■J-J Gaussian ensembles are used to model the chaotic part. A closed integral rep- 

?3 resentation for all correlations in the case of broken time reversal invariance is 

Ch derived by employing supersymmetry and the graded eigenvalue method. In 

particular, the two level correlation function is expressed as a double integral. 
For a correlation index, exact diffusion equations are derived for all three univer- 
O sality classes which describe the transition to the chaotic regime from arbitrary 

initial conditions. As an application, the transition from Poisson regularity to 
chaos is discussed. The two level correlation function becomes a double integral 
for all values of the transition parameter. 
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1 Introduction 



In classical and quantum mechanics, the interaction of a system is very often a 
sum of different contributions whose individual properties are well understood. The 
Hydrogen atom and the electron in a magnetic field, for example, are two of the 
best studied systems in quantum physics. The Hydrogen atom in a strong magnetic 
field, however, the combination of these two systems, could only fairly recently be 
investigated in full detail [1]. The total interaction is a sum of three competing 
terms, the Coulomb, the paramagnetic and the diamagnetic one. The infiuence of 
the different contributions can most instructively be seen when analyzing observ- 
ables that are sensitive to regular or chaotic behavior of the system. We use the 
designation "regular" or "chaotic" for a quantum system whose classical analogue 
is integrable and regular or non-integrable and chaotic, respectively. Observables of 
particular interest for those studies are the nearest neighbor spacing distribution and 
the spectral rigidity [2]. In the case of the Hydrogen atom in a strong magnetic field, 
these two observables approach the Wigner-Dyson predictions [2, 3, 4, 5] gradually 
as the magnetic field strength is increased, and thus the system makes a transition 
from a regular to a chaotic one [1]. According to this observation, the Hamilto- 
nian of a statistical model describing a regular system in an external field which 
induces chaotic motion can be written as the sum of a regular and a chaotic part. 
The theory of random matrices [3] provides the Gaussian Ensembles as well under- 
stood models for the chaotic part. Due to the general symmetry constraints, a time 
reversal invariant system with conserved or broken rotation invariance is modeled 
by the Gaussian Orthogonal (GOE) or Symplectic Ensemble (GSE), respectively, 
while the Gaussian Unitary Ensemble (GUE) describes the fiuctuation properties of 
a system under broken time reversal invariance. These ensembles provide a generic 
description of chaotic spectral fiuctuations. However, regular fiuctuation properties 
are not generic in this sense, they rather differ from system to system. One example 
is the harmonic oscillator whose energy levels are equidistant. Our main interest in 
this work will be the opposite situation, i.e. a regular system that does not show 
any spectral correlations. Such a system is then modeled by the Poisson Ensemble. 

Statistical models of this type can also be employed to understand similar trans- 
itions in other systems. Classical and quantum billiards [4] are very useful systems 
for chaos studies because of their relatively small number of degrees of freedom. 
Starting from a regularly shaped billiard, slight changes of the geometry often make 
the dynamics gradually or even abruptly chaotic. The Henon-Heiles system [4] 
can be viewed as the truncation of a regular and integrable system such that the 
transition from regular to chaotic behavior occurs as energy increases. Similar ob- 
servations have been made in molecular spectra. Hence, the decomposition into a 
regular and a chaotic part is not always dictated by externally tunable parameters 
like magnetic field or deformations of the geometry. Sometimes, a very complicated 
system becomes mathematically feasible if written as a sum of a somehow smooth. 
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dominating term and a correction. The nuclear shell model [6] might be among 
the most successful examples for the utility of this idea. The potential is a sum of 
an usually integrable mean field part and an often chaos producing residual inter- 
action. Another example from nuclear physics is provided by heavy ion reactions 
where compound nuclei are produced with very high angular momentum. Their 
rotational motion along the yrast line is collective and integrable. However, often 
there is an inner temperature due to single particle excitations creating chaos which 
becomes visible in a spreading of the electrical quadrupole transition strength [7]. 

Recently, it has been shown that the statistical random matrix model involving 
a sum of a regular and a chaotic part is even useful to analyze spectra in classical 
wave mechanics. The spectral fiuctuations of elastomechanical eigenmodes in Alu- 
minum blocks manufactured in the shape of three dimensional Sinai-billiards show 
a transition which closely parallels the transition from Poisson regularity to chaos 
in quantum mechanics [8]. 

Obviously, the detailed study of the transition from regularity, in particular Pois- 
son regularity, to chaos in such a random matrix model is a worthwhile task. Of 
the numerous studies addressing this problem, we mention particularly the work 
of French, Kota, Pandey and Tomsovic [9] and the work of Leyvraz and Selig- 
man [10]. In both studies, the two level correlation function for small values of the 
transition parameter is derived. Detailed numerical simulations can be found in 
Ref. [7]. However, despite several attempts, a full hedged analytical discussion of 
this transitions is still lacking. The calculation that comes closest to a solution of the 
mathematical problems is most likely the one presented by Lenz [11]. Surprisingly, 
other but related transitions like the breaking of time reversal invariance could be 
calculated analytically. In a statistical model, the Hamiltonian of a system with 
gradually broken time reversal invariance is expressible as a sum of two contribu- 
tions, one drawn from the GOE, the other from the GUE. The spectral fiuctuations 
of this transition were worked out analytically by Mehta and Pandey [12]. Using the 
level number variance as a measure, upper bounds for the time reversal invariance 
breaking in nuclear physics were given [13]. 

In this work, we study transitions from arbitrary to chaotic fiuctuation proper- 
ties. Our main interest is the transition from Poisson regularity to chaos, but many 
of our results are very general and include this transition and the time reversal in- 
variance breaking as special cases. The only general restriction we make is that 
the Hamiltonian can be written as the sum of two parts such that one of them can 
be modeled by a Gaussian Ensemble. Although most of the physically interesting 
transitions fall into this class, there are exceptions. The breaking of a quantum 
number, for example, has of course a completely different Hamiltonian structure. 
The infiuence of isospin breaking on the fiuctuation properties has been studied ex- 
perimentally [14] and theoretically [15] in great detail. Those systems are not the 
subject of this work. Moreover, for technical reasons, most of our results involving 
integral representations are derived for transitions toward the GUE. 
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For the mathematical discussion of the transition we rely on supersymmetry. In 
Refs. [16, 17] supersymmetric nonlinear ci-models were derived for a wide class of 
statistical quantum systems with chaotic dynamics. In Refs. [15, 18] a particularly 
suitable technique was developed in order to treat these models in the pure spectral 
case. 

In Section 2, we formulate the random matrix model and introduce our notation. 
Supersymmetric methods are applied in Section 3. In Section 4, we derive diffusion 
equations for the transition. The transition from Poisson regularity to chaos is 
worked out in Section 5. We discuss our results in Section 6. Certain calculations 
are performed in the Appendix. 

2 Transitions toward the GUE 

In Section 2.1, we summarize the main results for pure ensembles and introduce our 
notation. In Section 2.2 the model for the transition is introduced. For simplicity, 
all formulae are derived for the GUE. With the aid of Ref. [3], however, everything 
can easily be mapped onto the GOE or the GSE. 

2.1 Pure Ensembles 

We consider two ensembles of Hermitean N xN matrices H^'^^ and H^^\ respectively, 
with volume elements 

N 

d[H('^] = HdH^^^HdReH^^^dlmH^^^, j = 0, 1 (2.1) 

n=l n>m 

and normalized probability density functions Pj^\h^^'^)^ j = 0, 1. Some well known 
formulae which are true for arbitrary ensembles are given in the following. The 
A;-level correlation functions 

R^^\x„...,x,) = ^1 pO-)(^0-))-Qi^tr-^^^rf[i7W] (2.2) 

p=i P 

measure the probability of finding a level in each of the unit intervals around the en- 
ergies xi, . . . , Xk- The energies carry a small imaginary increment, x~ = Xp — ie, the 
limit e ^ has to be taken after the ensemble average. Notice that this definition 
differs slightly from the Dyson-Mehta definition [3] since it includes terms propor- 
tional to S(xp — Xq). We are not interested in those contributions and will therefore 
proceed as in Refs. [15, 18] by always assuming x,p ^ Xg for all p ^ q. In order to 
study the fiuctuations, the correlation functions have to be unfolded by measuring 
the energies on the scale of the mean level spacing D^^^ = l/i?|"'^(0), j = 0, 1 for in- 
finitely many levels. We define the dimensionless energies ^p = Xpj D^^\ p = 1, . . . ,k 



4 



which have to be held fixed in the limit N ^ oo. The unfolded A;-level correlation 
functions 

= ]im {D(^y)'^Ri'\x^,...,Xk) (2.3) 

iV ^oo 

measure the correlations on the energy scale with mean level spacing unity every- 
where. All physically important ensembles allow such a complete removal of the 
level densities from the higher correlations. The functions xj^''\^^ . . . ^ ^k) depend 
then only on the differences ^p — ^q. Those fiuctuation properties are called transla- 
tion invariant or generic. 

By including the real part of the Green function, we introduce the correlation 
functions 

R^^\x„...,x,) = ^1 pW(^0-))-Qtr-^^^rf[i7W] (2.4) 

p=i P 

where the imaginary increments can now lie on both sides of the real axis, we write 
x^ 



— ± ie. Due to the identity 



Imtr -T- = — (tr rr — tr— r ) , (2.5) 

xp - Hii) i2 \ xp - Hii) x+ - Hii) J ' ^ 

the physically relevant correlations (2.2) can always be constructed as a linear com- 
bination of the functions (2.4). As in Ref. [15], we use the operator symbol S in 
order to indicate this procedure. The functions (2.4) have the advantage of being 
expressible as the derivatives 



12.6) 



J=0 



of the generating function 



J p=i det {^xf - Jp- HO)j 

with source variables Jp, p = 1, . . .,k. For later purposes we write the argument as 
X + J with the 2k x 2k matrices 

X = diag(xi,xi,. . .,Xk,Xk) 

J = diag(-/i,+/i,...,-/fc,+/fc) . (2.8) 

The generating function is normalized, i.e. at / = we have zP(x) = 1. The 
function '^zj^^\x + J) generates the correlations (2.2). We emphasize again that all 
these formulae presented so far are true for arbitrary probability densities. 
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By choosing a Gaussian distribution for the matrices H^^\ 

with variance f ^, we make one of the ensembles a Gaussian Unitary Ensemble (GUE) 
with chaotic fluctuation properties. One usually fixes the energy scale by setting 
V2v = 1. We have the well known result [3] 

N-l 

KN(Xp,Xg) = J2 Vn(Xp)VniXq) (2.f0) 
n=0 

where Lpn{xp) is the ra-th oscillator wave function. The lower order correlations are 
removed in the A;-level cluster functions 

m 

Ti^\x,,...,x,) = J2i-ir-''im-iy. n 4i, , (2.11) 

G m'=l 

where G denotes a partition of the indices (1, . . ., A;) into subgroups (Gi, . . .,Gm)- 
The inverse reads 

m 

R^^Hx,,...,x,) = n 4i • (2.12) 

G m'=l 

Besides the unfolded A;-level correlation functions (2.3) we also define the unfolded 
A;-level cluster functions 

yi^)(ei,...,6) = lim (D('y)'^Ti'Hx^.....Xk). (2.13) 

iV ^oo 

The other probability density function pj^^H^^"^) we do not specify and thereby 
allow for arbitrary fiuctuation properties. 



2.2 Transition Ensemble 

We want to study the fiuctuations of the ensemble of matrices 

H(a) = + (2.14) 

as a function of the transition parameter a. In this formulation, the arbitrary en- 
semble is the initial condition, and as we increase a, the GUE fiuctuation properties 
will become stronger and eventually dominate. It is advantageous to absorb the 
transition parameter into the GUE by defining aH^^^ as new matrices. This can 
formally be achieved by replacing y^v in equation (2.9) with a such that a^/2 is 
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now the variance of the Gaussian distribution. We write pj^\H^^\ a) to indicate 
the dependence on the transition parameter. At the initial condition, the Gaussian 
distribution becomes a ^-function, 

lim pj^\H^^\ a) = 6(H^^y) . (2.15) 
The A;-level correlation functions of this transition ensemble are given by 



a) 

k 

Tm ir 

i?(o) - HW 



TTlmtr 1- (2.16) 

■p=\ -^p 



and the A;-level cluster functions Tk(xi, . . . , Xk, a) are defined as in Eq. (2.11). Again, 
easier to evaluate are the functions 

i?i(xi.....,n.o) = if /"<([if<°)]4°'(if<°))/"<([ifW]Pi"(5W.a) 



'p = l '-^ 



since they can be written as derivatives 



Id'' 



(2.18) 



j=o 



of the normalized generating function 



a 



n A ( • (2-19) 

pii det (4 -Jp- i?(o) - HW) 



Again, the physically relevant correlations (2.16) are generated by the function 
^Zk(x + J, a). Because of Eq. (2.15) we see immediately that the correlations of the 
arbitrary ensemble emerge for vanishing transition parameter, 

lira Rk(xi, . . .,Xk,a) = Rf\xi,...,Xk) 

lim Zkix + J,a) = zj^\x + J) . (2.20) 
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The limit a ^ oo has to yield the GUE result. In order to check that, we do not 
absorb a into H^^^ but use the Hamiltonian in the form (2.14), we find 

k ^ 

TTlmtr ■ — . (2.21) 

Xp/a- H(o)/a- H(^) ^ ' 

Obviously, the energies have to be rescaled when we compute the limit a ^ oo by 
holding fixed Xp = Xp/a. The factor in front of the integral takes care of the 

proper rescaling of the differentials, 

lim Rk(xi^ . . . ^ Xk^ a) dxi ■ ■ ■ dxk = R^i^\xi^ . . . ^Xk) dxi ■ ■ ■ dx^ (2.22) 

a^oo rev/ 

as required. Equivalently we can write 

lim a'' Rk(xi, . . .,Xk,a) = R^p (xi, . . . ,Xk) (2.23) 

a^OO /rev 

and similarly for the generating function. 

In the definition (2.14) of the transition ensemble, the parameter a can take 
values between zero and infinity. Since the ensemble is obviously symmetric in the 
transition parameter, we consider only non-negative values of a. In certain situ- 
ations, however, it is more convenient to choose alternative definitions. In general, 
we can write 

H{a) = 7(°)(a)i?(°) + j^^\a)H^^'> (2.24) 

where 7(°)(q!) should be symmetric and '•/^^'^(a) antisymmetric in a. A good example 
is 7(°)(q!) = Vl — ct^ and '•/^^'^(a) = a. Here, a varies only between zero and unity. 
One verifies easily that the limits a ^ and a ^ 1 give the correct results for the 
arbitrary ensemble and for the GUE, respectively. Another choice is 7(°)(q!) = cos a 
and 7(^)(q!) = sin a, where a varies between zero and 7r/2. We will mainly work 
with the definition (2.14). 

As in the case of the pure ensembles, we have to unfold the correlation functions 
for large level number N by removing the dependence on the level density. We meas- 
ure the fiuctuations on the scale of the mean level spacing D = 1/Ri(0,a), i.e. we 
map the spectrum onto an equivalent one with mean level spacing unity everywhere. 
Obviously, this can make sense only if the level densities of the pure ensembles in- 
volved cover the same section of the energy scale x. In other words, R^f\x) and 
R^^\x) should have the same asymptotic dependence on the level number. The GUE 
level density R^j^\x) becomes asymptotically the famous Wigner semicircle [3] which 
behaves as ViV"- Thus, the mean level spacing is proportional to l/ViV- We define 
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new energies = Xp/D, p = 1, . . . ,k. The transition parameter a is defined on the 
original scale and has therefore also to be unfolded. The new transition parameter 

A = a/D (2.25) 

was first introduced by Pandey [19], its universality was shown in many applica- 
tions [7, 9, 12, 15]. Since A has to be held fixed in the limit N oo, the original 
transition parameter a behaves formally as l/ViV, just like the mean level spacing. 
The A;-level correlation functions on the unfolded scale 

Xfc(ei,...,6,A) = lim D^Rk{xu...,Xk,a) (2.26) 

N^oo 

are translation invariant. By including the real parts of the Green function, we also 
define the functions Xk{S,, . . . , .^fc. A) for later purposes. 



3 Application of Supersymmetry 

In Section 3.1 the supersymmetric representation is introduced. The graded eigen- 
value method is applied in Section 3.2. In Section 3.3 a Pastur equation for the 
level density is derived. The unfolding procedure is performed in Section 3.4. In 
Section 3.5, the correlation functions on the unfolded energy scale are expanded in 
the transition parameter. Translation invariance and its consequences are discussed 
in Section 3.6. In Section 3.7, the level number variance is worked out. 



3.1 Supersymmetric Representation 

The average over the GUE in the generating function (2.19) can be evaluated by 
using the supersymmetric technique developed in Refs. [16, 17]. The advantage of 
this technique is a dramatic reduction of the degrees of freedom which becomes 
obvious in the decoupling of the number of integrations to be performed from the 
level number N . The procedure can straightforwardly be applied to our case and 
will therefore not be discussed here. Since we will later not employ the coset method 
of Refs. [16, 17], the reader who is not familiar with supersymmetry in this context 
is advised to study Ref. [18] which also defines the notation for the remainder. 

The GUE is mapped onto an ensemble of 2k x 2k Hermitean graded or super- 
matrices a. In pg-block notation [17, 18], these matrices are defined as kxk matrices 
whose entries apg are 2x2 supermatrices of the form 



PQ PQ 



f3.1) 
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where (J^q^c =1,2 and are commuting and anticommuting variables, respect- 
ively. The Cartesean volume element reads 

2 k 

d[a] = n n n dRea^^q^dlvaa^^q^ H ^« • (3-2) 

c=l p=l p>q p,q 

There is a normalized graded probability density function 

Qk{a,a) = 2'=('=-i) exp (^-^trgd^^ (3.3) 
having the well defined limit 

l\mQkia,a) = S{a) (3.4) 

where the right hand side is a properly defined ^-function in superspace [16, 18]. 
The generating function (2.19) can then be expressed in the form 

Zk{x + J,a) = J Pj^\H(^y) J d[a] Qk{<T, a) 

detg-i ((a;± + / - d) ® Ijv - lafc ® (3.5) 

in which we introduced the N x N and the 2k x 2k unit matrices Ijv and l2k- Notice 
that the order k of the correlation is now directly refiected in the dimension 2k 
of the supermatrix a, this causes the enormous simplification. For the two-level 
correlation function of the Gaussian Orthogonal Ensemble (GOE), a similar integral 
representation of the generating function was obtained in Ref. [7]. 

3.2 Graded Eigenvalue Method 

In order to further evaluate the generating function we use the graded eigenvalue 
method which was developed in Ref. [18]. After a simple shift of variables and a 
interchange of the order of integration we find 

Zk(x + /,«) = J d[a] Qk(a^ - x - .J,a) j d[H^^\ Pj^HH^''^) 

detg-^ ((7®ljv-l2fc®i?(°)) (3.6) 

such that due to the direct product structure of the argument the superdeterminant 
depends only on the eigenvalues of a and H^'^K All angular degrees of freedom are 
now in the graded probability density function. We therefore go over to eigenvalue 
and angle coordinates 

a = u~^su (3.7) 
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where u is an element of the unitary supergroup U(k/k). The diagonal matrix of 
the eigenvalues 

s = diag(sii, isi2, . . . , Sfci, «Sfc2) (3.8) 

contains k eigenvalues Spi of boson boson and as many eigenvalues isp2 of fermion 
fermion type [17, 18]. 

For the integration over the supermatrix, we need the transformation of the 
Cartesean volume element (3.2) to eigenvalue- angle coordinates [18], 

d[a] = Blis)d[s]di2iu) , (3.9) 

where d[s] is the product of the differentials of the eigenvalues and djiiu) is the 
invariant Haar measure of the unitary supergroup. The square root of the Jacobian, 
here referred to as Berezinian, is given by 



Bk{s) = ^zj— : = det 



1 



Spl - lSq2 



(3.10) 

p,q = l,...,k 



which reflects [18] the determinant structure (2.10) of the GUE correlation functions. 

The angular average over the unitary supergroup can be performed by using 
the supersymmetric generalization [18] of the Harish-Chandra-Itzykson-Zuber in- 
tegral [20, 21]. We define a second Hermitean supermatrix v~^rv of the type (3.7) 
where v is superunitary and r diagonal. The angular average over the shifted Gaus- 
sian probability density function yields 

Qk{u su-v rv,a)dfi{u) = ^ , , . (3.11) 

Bk[s)-tSk[r) 

where the Gaussian kernel on the right hand side is given by 

Gk{s,r,a) = 'iet[g{spi - rqi,a)l^^g^-^^_^^ det [^((5^2 - ?-g2, a)]p^g=i^...^fc 

g(z,a) = ^_ exp ( ) . (3.12) 

This result does not depend on the superunitary matrix v because of the invariance 
of the Haar measure djiiu). Since the determinants and the Berezinian (3.10) are 
antisymmetric, simple interchanges of variables show that all terms in the Laplace 
expansion of Gk{s, r, a) contribute the same if further integration over the eigenval- 
ues s or r is required and the rest of the integrand has the same symmetries. Hence, 
in our case, it is sufficient to write 

Gk{s,r,a) = Gk{s-r,a) = exp f — ^trg(s - r)^ j (3.13) 
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under the integral. Furthermore, the integration over s requires to take a new type of 
boundary contributions [16, 17, 22, 23] into account which do not occur in ordinary 
analysis. Following Ref. [18] we have to add the term (1 — r](r))S(s)/ Bl(s) to the 
right hand side of the formula (3.11) which contains the distribution 



In our particular case, we want to apply formula (3.11) to the integral (3.6). Thus, 
we have a diagonal shift matrix, i.e. f = 0, with eigenvalues given by r^i = Xp±ie — Jp 
and irp2 = Xp ± ie + Jp. It is then obvious from Refs. [18, 24] that formula (3.11) 
is valid for all choices of the signs of the imaginary increments. Consequently, the 
procedure to go over to the correlations of the imaginary parts which was designated 
by the operator symbol S does still work. 

After performing the angular integration we shift the imaginary increments back 
into the superdeterminant. A comparison with the definition (2.7) shows that the use 
of the eigenvalue angle coordinates allows to express the average over the arbitrary 
ensemble as nothing else but the generating function of the arbitrary correlations. 



for the generating function of the transition ensemble. The boundary contribution 
ensures the normalization at / = such that Zk(x,a) = 1. However, in the limit of 
vanishing transition parameter a, the Gaussian kernel (3.12) becomes a ^-function, 
in particular we find the diagonal term S(s — x — J) for the expression (3.13). The 
integral yields then the initial condition Z^\x + /) for all values of the source 
variables, including / = in perfect agreement with Eq. (2.20). Hence, in this 
special case, the boundary contribution 1 — ri{x + /) is obsolete. We keep that in 
mind, but we do not indicate that through a special notation. There is still another 
caveat in order. It was shown in Ref. [25] that the treatment of the boundary terms 
which is used here can miss some contributions to the real parts of the correlation 
functions. Nevertheless, since everything becomes correct in the physically relevant 
limit of infinitely many levels and since we are only interested in the imaginary parts 
anyway, we can proceed with formula (3.16). 

It is now very easy to go over from the generating function (3.16) to the A;-level 
correlation functions. For a = we find by construction 




(3.14) 



Zf\s) = I Pj^HH(^y) detg-1 [s^ ® Ijv - hk ® d[H(^)] (3.15) 



where s takes the place of a; + /. Collecting everything we arrive at 
Zk{x + J,a) = 1 - ri{x + J) + 




j Gk{s -x-.J,a) Zf\s) Bk{s)d[s] (3.16) 



a;fc,0) = Rf\xi,... 




(3.17) 
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which is consistent with Eq. (2.20). For non-zero a, we arrive at 



-1 



TT 



Rk(xi, . . .,Xk,a) = '—J- I Gk(s - x,a) zj^^s) Bk(s)d[s] (3.18 



by using the methods of Ref . [18]. In order to find the correlations of the imaginary 
parts we simply have to replace zj^\s) in the integrand of Eqs. (3.16) and (3.18) with 
?yZ|,°^(s). The A;-level correlation function for the transition from arbitrary initial 
conditions toward the GUE is represented as a 2A;-fold integral. Formula (3.11) 
is crucial for the graded eigenvalue method. In a single step, all angular degrees 
of freedom including all Grassmann variables are integrated out. This is true for 
arbitrary index k of the correlations. In the coset method used in Refs. [16, 17], the 
evaluation of correlations with indices k > 2 becomes increasingly complicated. 

3.3 Level Density and Pastur Equation 

The graded eigenvalue method is especially suited for a discussion of the correlations 
in the case of transitions [15]. However, the asymptotic behavior of the level density 
for large level numbers can easier be found by using the saddlepoint method of 
Refs. [16, 17]. Starting from Eq. (3.6) and using Eq. (3.15) we can write 

Zi(x + J,a) = J exp L(a, X + J, a) d[a] (3.19) 

with the Lagrange function 

L(a,x + J,a) = -^trgia - X - Jf + lnZ[°\s) . (3.20) 

Since at the saddlepoint we can set / = 0, the saddlepoint value a = ao'is determined 
through the equation 



= dL(a,x,a] 



trg 



^{a-x)+^lnZf\s) 
a'^ da 



da (3.21) 

J (7 = (7o 



which is solved in Appendix A. The solution is diagonal, both eigenvalues are equal, 
= is2o, and satisfy the implicit equation 
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^(sio-a;) + 7r4°)(sio) = (A ^ ^) (3.22) 

where we omit the index A; = 1 in the eigenvalues and the energy. Here, of course, 
-Ki°^(sio) stands for the leading term in an asymptotic expansion for large level 
number A. 

The asymptotic form of the one point function is then found to be 

Ri(x^a) = (x — sio(x^ aj) = R^^\sio(x^ aj) (3.23) 
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where sio = sio(a;,Q!) is the solution of the implicit equation (3.22), details are 
sketched in Appendix A. The level density 

Ri(x^a) = -lmsio(x^a) = Imi?i°Vsio(a;, a)) (3.24) 

is the imaginary part of the one point function. All these results have to be under- 
stood asymptotically, i.e. for large level number N. The level density Ri(x,a) is 
always normalized if the arbitrary initial condition is normalized to the total level 
number, 

+00 +00 

Ri(x^a)dx = / R^f\sio(x^ aj) dx 



-00 —00 



+00 ^ +00 



/(O)/ X , ■'TO! f ^(0)/ X dR] 

R\ '{siojdsio H ^Im J R\ '{sio)— — (sio) «sio 



-00 



= N . (3.25) 

We changed the variable of integration from x to sio with the help of the Eqs. (3.22) 
and (3.23). The second integral in the second line is zero since for consistency 
reasons the level density R^f\x) has to vanish for infinitely large energy argument. 

Pastur [26] derived the implicit equation (3.22) without using supersymmetry. It 
appears here through a simple saddlepoint approximation because supersymmetry 
decouples the level number N from the integration. The implicit equation for the 
pure Gaussian ensembles can easily be obtained by choosing the unaveraged Green 
function 

Rf'\sio) = -— (3.26) 

TT Sio 

as initial condition and by setting a = 1, we find from Eq. (3.22) 

2isw-x) + — = (3.27) 

•SlO 

which yields the well-known semicircle law. More interesting, however, is the case 
of having a semicircle as initial condition. 



i?;°)(sio) = i (^sio + i^2A - s2^j . (3.28) 
We solve equation (3.22) for sio and find from formula (3.23) another semicircle law 



Ri(x,a) = — , ^ I , ^ + i\ 2N - ( , ^ 1 I (3.29) 



which involves the rescaled energy x/\/l + . Note that this result is correct for 
all values of a. 
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3.4 Unfolding Procedure 

The main advantage due to the application of supersymmetry and the graded eigen- 
value method is the decoupling of the degrees of freedom of the supermatrices from 
the level number which appears only in the initial condition zj^\s). This feature 
makes it possible in many situations to perform the unfolding procedure in the gen- 
erating functions. Besides the new energies = Xp/D we also define new source 
variables jp = Jp/D ,p= 1, . . .,k and the 2k x 2k matrices ^ and j analogously to 
Eq. (2.8). The A;-level correlation functions on the unfolded energy scale depending 
on the transition parameter A can then be written as the A;-fold derivative 



Id'' 



(3.30) 



of the unfolded generating function 



Zk(i + j,>^) = lim Zk(x + J,a). (3.31) 

It is useful to unfold the integration variables s in Eq. (3.16) by making the change 
of variables s s/D. We find 

^fc(e + i,A) = i-vi^+j) + 

^ ^ Gk{s-C-J,>^)z'f\s)Bk{s)d[s] (3.32) 



where the unfolded initial condition is given by 

4°)(.) = lim zj^HDs) . (3.33) 

iV ^oo 

Consequently, the whole unfolding procedure is naturally shifted into the initial 
condition. As we will see in the following two sections, this is very useful from a 
conceptual viewpoint and also for applications. Analogous to Eqs. (3.17) and (3.18), 
we find for the A;-level correlation functions on the unfolded scale 

Xfc(ei,...,6,0) = Xf)(ei,...,6) (3.34) 
for A = and by using Eq. (3.30) we have 

Xfc(ei,...,6,A) = J Gk{s-^,\)z';^\s)Bk{s)d[s] (3.35) 

for non-zero A. As on the original scale, we simply have to replace z^i^\s) by ?y2;|,°^(s) 
in order to go over to the physically interesting contributions from the imaginary 
parts of the Green functions. 
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A comment on the alternative definitions (2.24) of the transition parameter a 
is in order. The parameterization 'y^'^'^(a) = \J\ — a"^ and 7(^^(0!) = a gives on 
the unfolded scale the same result (3.32) since a behaves formally as l/-\/iV" and 
thus limjv^oo 7^°H«) = 1- This is also true for the choice 'y^'^'^(a) = cos a and 
7(^)(q!) = sin a. We conclude that a wide class of parameterizations yields the same 
result (3.32) on the unfolded scale. The limit A ^ gives the fiuctuations of the 
arbitrary initial condition. The GUE results have to show up for A ^ 00. However, 
it is known from several applications [7, 9, 12, 15] that the chaotic features start 
dominating the fiuctuations for much smaller values of A. In most cases, a value of 
A around unity will be enough to spread the GUE correlations over the spectrum. 

3.5 Expansion in the Transition Parameter 

The A;-level correlation functions Xk(^i, . . . ,^k,X) following from Eq. (3.35) have a 
remarkable scaling property which shows up when the variables s/A are introduced 
for the integration, 

Xfc(ei,...,6,A) = ^ J Gk{s-a\A) '^'\!:^'^ Bk{s)d[s] . (3.36) 

The transition parameter A has disappeared from the variance of the Gaussian func- 
tion and occurs only in the initial condition and in the combination ,^/A. Since this 
rescaling of the energies is uncritical for all parameter values, we may hold the ratios 
^p/X,p = 1, . . .,k fixed and expand the initial condition in the transition parameter. 
Two boundary conditions are present: firstly, the correlations must obviously be an 
even function in A and secondly, there should not be a singularity for vanishing A. 
Hence, it is reasonable to assume an expansion of the form 

^^k (AgJ _ (0). X 

xk - 2^ 

for the initial condition which yields the expansion 

Xfc(ei,...,6,A) = ^ — Xfc,„(ei/A,...,6/A) (3.38) 

for the A;-level correlation functions. The coefficient functions 

Xfc,„(ei/A,...,6/A) = jGk{s-a)^,l)vj'^^l{s)Bk{s)d[s] (3.39) 

depend only on the rescaled energies ^p/X. The expansion (3.38) generalizes an 
observation made by Leyvraz and Seligman [10]. For a Poisson initial condition, 
they calculated the two level correlations in zeroth order by using a self-consistent 
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(3.37) 



perturbation theory and found that it depends only on (,^1 — ^2)/^- Formula (3.38) 
shows that the coefficient functions for all indices k of the correlation have such a 
scaling property provided the form of the expansion (3.37) is valid. An important 
caveat is in order: If the arbitrary ensemble has correlations, the zeroth order term 
of the expansion has to reproduce them. Since they are not on the scale ^/X, 
there should be an additional term to the zeroth order. However, this term will be 
specific for the ensemble under consideration and can therefore not be discussed in 
general. With the same caveat, the same scaling property can easily be derived for 
transitions toward the Gaussian Orthogonal or Symplectic Ensemble, GOE or GSE, 
respectively, again generalizing the corresponding observation made by Leyvraz and 
Seligman. 

3.6 Translation Invariance and its Consequences 

The correlation functions (3.35) are translation invariant as required provided the 
initial condition z^i^\s) depends only on the differences of the eigenvalues s. In order 
to show this, we shift the integration variables s by the matrix ^ in the integral (3.35). 
If z^i^\s) depends only on the differences of the eigenvalues, the correlations are 
obviously only functions of the differences — ^q. This is exactly the translation 

invariance. Of course, the generating function z^i^\s) will depend on the differences 

of the eigenvalues if the the arbitrary correlations xj^\^^ . . . ^ ^k) are translation 
invariant. Thus, the correlations of the transition are translation invariant if the 
arbitrary initial condition is translation invariant. 

The fact that translation invariance is visible in the integral representation (3.35) 
has some helpful consequences. For the level density, it is shown in Appendix C that 
Eq. (3.35) directly implies 

Xi(e,A) = 1 (3.40) 

for all values of A. If we start from a properly unfolded arbitrary spectrum, the level 
density will always be unity during the entire transition, as it should be. 

In the case of the two level correlations, translation invariance allows to reduce 
the four fold integral (3.35) to a double integral. We introduce the sum r = + ^2 
and the difference r = ,^1 — ,^2 of the energies. According to this rotation in energy 
space, we make the change 

si = sii + S21 , h = sii - S21 , 

S2 = S12 + S22, t2 = S12 - S22 (3-41) 

of the integration variables. Since z^2^\s) depends only on the differences of the 

eigenvalues, it cannot depend on si + is2. We write Z2^\s) = Z2^\si — is2,ti^t2) 
which, for reasons of consistency, should be even in each of the differences ti and ^2- 
As shown in Appendix D, the integrals over si and S2 can then be performed and 
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the dependence on r disappears. The two level correlations can thus be cast into 
the form 



+00+00 



Ut 



11-2 



r.2^,2^2^4°V,^l,^2)rf^lrf^2 (3.42) 
I'l + '■2) 

which is usually not amenable to further analytical treatment. For the higher cor- 
relations with k > 2, similar simplifications are likely to exist. 

The expansion in the transition parameter discussed in the previous subsection 
can be written more explicitly for the two level correlations using the translation 
invariant form (3.42). Introducing ti/X and t2/X as new integration variables we 
find 

+00+00 

X2(r,A) = J exp(^-^{tl + tl)^ smhjhsmjt2 

— 00 — 00 

ht2 s4°^(0,Aii,Ai2) 



A2 



dtidt2 . (3.43) 



We hold the rescaled energy difference r/A fixed and expand the initial condition in 
the transition parameter, 



s4°V,Aii,Ai2) 



E 



X2n 9'^+lcj^(0)(0,Aii,Ai2) 



A2 fo(^+^y- ^(^')'^+' 

which gives inserted into Eq. (3.43) the expansion 

00 -1 2ra 

^2(r,A) = ^— X2,.(r/A) 



(3.44) 



A2=0 



n=0 



(3.45) 



for the two level correlations. The coefficient functions 

+00+00 

exp 



^2,n(r/A) 



TT^in + 1) 



1 



2^^! + ^2) ) sinh -ti sin -^2 



tit2 9'^+^s4°V,Aii,Ai2) 



Hi + 4 



9(A2)'^+i 



dtidt2 (3.46) 



A2=0 



depend only on one variable, the rescaled energy difference r/A. Of course, the 
caveat at the end of the last subsection is also valid here. 
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3.7 Level Number Variance 



The variance of the level number on the unfolded energy scale studied as a func- 
tion of the length L of the energy interval under consideration is of great interest 
for the analysis of data since it gives information about the long range spectral cor- 
relation. It describes two level correlations and can therefore be expressed involving 
solely the two level cluster function [2]. This is also true for transitions and we can 
thus make use of the formula 

L L 
S2(i,A) = L - 2j(L-r)Y2(r,X)dr = -L(L-l) + 2 j (L-r) X2(r, X) dr (3.47) 



where we employed the relation X2(r, A) = 1 — Y2(r, A) between the two level correl- 
ation and cluster function. Using Eq. (3.42) we can write down an integral repres- 
entation for the level number variance for the transition from arbitrary correlations 
to GUE correlations. Since the r integration can be carried out trivially, we arrive 
at the double integral 

2 -l-oo-l-oo 

^\L,X) = -L{L-l)+'-^l I e^J-^(tl + tP^ '''' 



— (yj — cyj 

f . ,Lh . Lt2 , 2ht2 / , Lh Lt2 

smh -— sm -— -| — ^ tt 1 — cosh cos 



V A2 A2 tl + tl V A2 A2 



^-2 

for all values of the transition parameter. 



'^z^°\0,ti,t2)dtidt2 (3.48) 



4 Transition as Diffusion 

In Section 4.1, a general diffusion equation for the GUE is derived. The generaliz- 
ation to the GOE and the GSE is performed in Section 4.2. From these diffusive 
processes, stationary equations for the pure ensembles are constructed in Section 
4.3. In Section 4.4 these equations are formulated on the unfolded energy scale. 



4.1 Diffusion Equation 

In the Eqs. (3.6) and (3.16), the generating function Zk(x + J, a) is expressed as an 
integral over the initial condition zj^\s) and a Gaussian kernel in a Cartesean and in 
a curved space, respectively. This representation of the generating function allows 
to interpret the transition from arbitrary toward GUE correlations as a diffusive 
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process. The squared transition parameter a^, i.e. the variance of the Gaussian 
distribution provides a fictitious time 

t = (4.1) 

that controls the diffusion. The combination a; + / of the energies and the source 
variables plays the role of fictitious coordinates 

r = X + J (4.2) 

that define the space in which the diffusion takes place. We view these coordinates 
as the eigenvalues of a Hermitean supermatrix cifc which is diagonalized by a unitary 
supermatrix such that cifc = u^^^rufc- There is a gradient d/da^c and a Cartesean 
Laplace operator 

A = trg^ (4.3) 

in this space [18, 35, 36]. The graded probability density function (3.3) written in 
the form 

gfc(cTfc,0 = 2'=('=-i) exp (^-itrg4^ (4.4) 

satisfies the definition of a diffusion kernel, 
1 d 

-AQk(aic,t) = -Q^Qk((Tic,t) and limQk(aic,t) = ^(o-fc) . (4.5) 
Therefore, the integral (3.6) is equivalent to the diffusion equation 

-AZkir,t) = —Zkir,t) with lim Zfc(r, i) = • (4-6) 

Since Zk(r,t) depends only on the eigenvalues r, the diffusion takes only place in 
this curved space and the radial part [18] of the Laplace operator is sufficient. It 
can be written as 

A, = ——-^_^.Blir)i-_, (4.7) 

Bl(r) dr " ' dr ' ' 

where we introduced the vector r = (rn, • • •,rfci,ri2, • • •,'''^2) and the corresponding 
gradient d/drfor notational purposes. Notice that we choose the metric such that 
there are no imaginary units in front of the variables rp2 and that the dot indicates 
the standard scalar product. Hence, we arrive at the diffusive process 



1 d 

-ArZk(r,t) = —Zk(r,t) with limZfc(r,i) = 



for the generating function of the transition correlations where the generating func- 
tion of the arbitrary correlations is the initial condition. The level number N does 
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not appear in the diffusion equation itself, it is only implicitly contained in the initial 
condition. This means that asymptotic solutions for the generating function of the 
transition correlations can be found if the large N behavior of the initial condition 
is known. 

The diffusion is a physically very instructive picture. A Gaussian kernel propag- 
ates arbitrary correlations diffusively toward chaotic correlations. This formula- 
tion is completely equivalent to both integral representations, the original Carte- 
sean (3.6) and the one in the curved eigenvalue space (3.16). The knowledge of the 
formula (3.11) is not necessary for the derivation of the diffusion equation (4.8). 
In fact, formula (3.11) was derived, in the ordinary [21] and in the supersymmet- 
ric case [18], by using a diffusion equation of the type (4.8). There, however, the 
fictitious time had no physical meaning, here, it is the transition parameter. The 
fictitious coordinates (4.2) combine energies and source variables, the relevant space 
is therefore not just given by the k energies. The k source variables are necessary 
for a construction of all together 2k independent variables r which is essential for 
the procedure. The correlations functions (2.16), however, depend only on the en- 
ergies. Hence, there is no diffusion of the type (4.8) for the correlation functions 
themselves, only for their generating functions. The existence of the diffusion pro- 
cess in the space of the fictitious coordinates (4.2) proves once more how well suited 
supersymmetry is in the theory of random matrices. 

The special form of the Laplacean (4.7) allows the reformulation [18] of the 
diffusion in the curved space in a Cartesean space. A direct calculation shows that 
the function 



Wk(r,t) = Bk(r)Zk(r,t) with lF^"\r) = Bk(r)Zl^>(r) (4.9) 



which could serve as a starting point for a power series expansion in the fictitious 
time t if the integral representation (3.16) was unknown. 

A comment on the boundary contribution in formula (3.16) is in order. One 
might question that the distribution 1 — r](r) satisfies Eq. (4.8). However, since 
Zk(r,t) in the form (3.6) and the integral in the curved space representation (3.16) 



as the initial condition solves the diffusion problem 



The processes (4.8) and (4.10) have the formal solutions 





(4.11) 
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solve the diffusion equation (4.8), 1 — r](r) must solve it as well. This can also be 
shown in a direct calculation using the integral representation of r](r) which was 
constructed in Ref. [f8], one easily finds Ari](r) = as it should be since r](r) is not 
time dependent. In this sense, the boundary contribution I — r](r) which ensures the 
normalization of the generating function is a stationary point of the diffusive process. 
Nevertheless, a minor problem arises when the transition parameter vanishes. As 
discussed already in Section 3.2, the boundary contribution is obsolete at i = 0. A 
careful investigation of the power series in t resulting from the formal solution (4.11) 
proves that this new type of singularity is an artifact caused by the incomplete 
treatment of the boundary contributions of the real parts of the Green functions. 
We therefore do not worry about this further since we are only interested in the 
imaginary parts and the limit of infinitely many levels which are both known to 
come out correct [25]. It should be stressed that Rothstein [22] showed a cleaner 
way of treating the boundary contributions by modifying the integration measure 
in the curved space such that it still has the symmetries of the original, in our 
case Cartesean, space. Unfortunately, so far it was only possible to apply this 
technique [25] to the one point function, i.e. to 2 x 2 Hermitean supermatrices. 

4.2 Orthogonal and Symplectic Case 

Besides the initial condition, the only input for the derivation of the diffusion equa- 
tion is the observation that the Gaussian probability density (2.9) is mapped onto 
a graded probability density (3.3) which is as well Gaussian. In Cartesean coordin- 
ates, this feature is also present [16, 17] in the Gaussian Orthogonal (GOE) and in 
the Gaussian Symplectic Ensemble (GSE). Hence there exists a Cartesean integral 
representation of the type (3.6) for the generating functions 

Zfjkix + J,t) = j Qf^^{cj-x-.J,t)zfi{s)d[(j] (4.12) 

of the correlations in the case of transitions toward all three of the Gaussian En- 
sembles. The parameter [i takes the values 1, 2 or 4 for the orthogonal, unitary or 
symplectic symmetry class, respectively. Here and in the following two subsections, 
all probability density functions and all correlation functions and their generating 
functions have the additional index [i which indicates a proper adjustment of the 
definitions in Section 2 to the different symmetry classes. The variance of the Gaus- 
sian distribution for the three values of (3 is often written as 2f ^//3 with Eq. (2.9) as 
special case for /3 = 2. Thus, fixing the energy scale through the condition \/2v = 1 
yields an effective variance of 1//3. For all three symmetry classes we use the defin- 
ition (4.1) of the fictitious time t. Hence, the variance of the graded probability 
density function reads 2t/(3. For the GOE and the GSE, the supermatrices a have 
to have dimension 4k x 4k in order to implement the additional symmetry con- 
straints. The diagonalization analogous to (3.7) involves properly defined 4k x 4k 
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matrices u and the eigenvalue matrices 

s = diag(siii,si2i,isi2,«si2, • • •,Sfcii,Sfc2i,«Sfc2,«Sfc2) (GOE) 

s = diag(sii,sii,isii2,isi22, • • •,Sfci,Sfci,«Sfci2,«Sfc22) (GSE) (4.13) 

with 3k independent variables. However, the 4kx4k diagonal matrices of the energies 
and the source variables 

X = diag(a;i, si, si, si, • • • , a;/;, a;/;, a;/;, a;/;) 

/ = diag(-Ji,-Ji,+Ji,+Ji,- ■ ■ ,-Jk,-Jk,+Jk,+Jk) (4.14) 

contain only k independent variables each. The combination x + J has only 2k 
independent variables and cannot serve as the fictitious coordinate space. We have 
to make the replacement 

X + J — > r (4.15) 

where r has the same form as s in Eq. (4.13). For notational purposes, we introduce 
the vectors 

r = (riii,ri2i,ri2, • • •,rHi,rfc2i,rfc2) (GOE) 

r = (rii,rii2,ri22, • • •,'rfci,'rfci2,'rfc22) (GSE) (4.16) 

and the corresponding gradients. Moreover, by writing cifc = u^^^rufc with cifc and 
Mfc of the same form as a and u, respectively, we may define a Cartesean Laplace 
operator analogously to Eq. (4.3). Obviously, the diffusion takes only place in 
the curved space of the eigenvalues r and the radial part 

1 ^ O 

Bfjkir) dr ' ' ' dr 

suffices. The Berezinians of the transformation to eigenvalue angle coordinates can 
be worked out as in Ref. [18], 

„ t . np=ij>/ kpji - vii np>gj,/ kpji - ^/ii np>g(^?-p2 - ir^2f 
''(""^ " fr~f7T^^777p 

= n,„„(^„--„.F 

for the GOE and the GSE and B2h{r) = Bl(r) for the GUE. Thus, collecting 
everything we arrive at the diffusion equation 

^AfjrZfjkir,t) = j^Zpk{r,t) with lira Zpk{r,t) = zJfcV) (4-19) 
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for all three symmetry classes and for arbitrary initial conditions. In case of the GOE 
and the GSE, the fictitious coordinates r define according to Eq. (4.15) an extended 
space in which the generating function is propagated. In order to calculate the 
correlation functions, we have to restrict this space to the subspace defined by the 
combination x + J. The formal solution of the diffusion equation 



Zfjkir,t) = exp(^Afjr) zf^,{r) (4.20) 



V4 

allows an expansion in powers of the fictitious time t for all three values of /3. In 
order to calculate the correlation functions up to ra-th order in i.e. in a^, one has 
to evaluate AJ^[,Z^°^(r), ra' = 0, . . . , ra at a; + /. 

It should be emphasized that the diffusion equation (4.19) was derived without 
performing the angular average, i.e. without knowing the explicit solution of the 
supersymmetric Harish- Chandra- Itzykson-Zuber integral 

Tpk{s,r,t) = Qi3k{u~'^su - v~'^rv,t)dfi{u) (4.21) 



in the cases /3 = 1 and (3 = 4. Unfortunately, the explicit solution is only known in 
the unitary case (3.11). Even for the ordinary orthogonal and symplectic groups, this 
integral has not been evaluated yet. In other words, we cannot give an eigenvalue 
integral representation of the type (3.16) for transitions toward the GOE and the 
GSE. So far, the solution (4.20) is the only available one. This is closely related 
to the structure of the radial part (4.7) of the Laplacean. For /3 = 1 and /3 = 4, 
there is no ansatz in the spirit of Eq. (4.9) that would reduce the diffusion in the 
curved space (4.19) to a diffusion in a Cartesean space analogously to Eq. (4.10). As 
mentioned already, the explicit solution (3.11) of the group integral was derived [18] 
by using this property of the radial part in the unitary case. 

However, for the sake of completeness, we summarize the general formulae in- 
volving the kernel (4.21). By construction, it has to satisfy the diffusion equation 

^AprTpk(s,r,t) = ^^Tpk{s,r,t) (4.22) 
with the initial condition 

lirn r^fc(s, r, i) = j S(u~^ su — v~^rv) dfi(u) 

det[S(spi - r^i)]p,g=i,...,fcdet[^(gp2 - ?-g2)]p,g=i,..,fc 

and, moreover, it has the property 

Tfjkir,r',t + t') = [Tfjkir,s,t)Tfjkis,r',t')Bfjkis)d[s] . (4.24) 
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Thus, we can write the solution of Eq. (4.19) in the form 

Zfjkir,t) = j Vpk{s,r,t)zfi{s)Bpk{s)d[s] (4.25) 
which of course coincides with Eq. (3.16) in the case /3 = 2. 
4.3 Stationary Equations for the Pure Ensembles 

The diffusion equation (4.19) allows for arbitrary initial conditions. As a special 
case, we can also study the diffusion of the generating function of a symmetry class 
into itself. More precise, for a given /3, we take the generating function of the 
correlations of this symmetry class as the initial condition of the Eq. (4.19), 

= (4.26) 

where, as in Section 2.2, the upper index (1) indicates the pure Gaussian Ensemble. 
This situation is somewhat stationary since the propagation changes neither the 
symmetry nor the structure of the correlation. Of course, the energy scale is affected 
and becomes a function of the fictitious time. 

It is very easy to determine this rescaling from the definitions (2.16), (2.17) 
and (2.19) of the correlation and the generating function without using the diffusion 
equation or any explicit form of the correlation functions. Both probability density 
functions Ppj^(H^^^) j = 0, 1 are now Gaussian. Again, we move the transition 

parameter a into the standard deviation of P^]^(i?(i)). By defining = + 

H^^^ as new matrices for the integration, the H^'^^ integral becomes the convolution 
of the two Gaussians with variances 2f^//3 and 2v'^a'^/(3 which gives a Gaussian 
with variance 2f^(l + a^)/ (3. By using the elements of the matrices H^^^ j^Jl + 
as new integration variables, we see that the energies and source variables have to 
be rescaled by Vl + ofl ■ As always, we set ^/2v = 1. Thus, in a given symmetry 
class, we find the scaling law 

%.(.!,...,..,«) = ^==^<^(7=^'---'7l=f) (4-27) 



for the transition of the correlations into themselves. It follows just from the defin- 
ition of the correlation functions. The knowledge of the explicit solutions [3] for 
the pure ensembles is not necessary to derive Eq. (4.27). Also only employing the 
definition, we arrive at the scaling law 

,yrT2? 



Zpkir.t) = Z^^>{^==) (4.28) 



for the generating functions. Here, we use the fictitious time t = a^/2 and the 
fictitious coordinates r. 
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We now take advantage of the fact that the generating function satisfies the 
diffusion equation (4.f9). For our special choice of the initial conditions, we have 



with the trivially fulfilled initial condition 



lini4iM^=^ = Zj^ir). (4.30) 



We set ri = rl\J\ + It and = rl\J\ + It and rewrite the time derivative 

dt ^''^ ~ dn ' dt ~ i + 2t' dn ■ ^ ^ ^ 

Since the Laplacean obeys A^^ = ^/3rt/(l + 2^) we find 

-^^f^nZfiin) = -n.^^zfiin) (4.32) 

which contains no explicit time dependence anymore. Hence, the generating function 
satisfies a stationary partial differential equation in the rescaled fictitious coordinates 
Tf. Therefore, we drop the index t and 

^A,. + f.^.) Zfiir) = (4.33) 

is our final result. The boundary condition is of course the normalization zj^j^(x) = 
1. This stationary equation is equivalent to the original integral definition of the 
type (2.7) for the pure ensembles. Notice that the only input is the definition of the 
correlation functions and the existence of the diffusion equation (4.19) for arbitrary 
initial conditions. It is easy to rewrite the stationary equation as 



(1), 



r 



4 5,-.+ n-^/^^(0^^ = (4.34) 



which has an integrating factor such that 



■ exp [+^tTgr'j B^k(r)^^ = . (4.35) 

From this equation, a new integral representation for zj^j^(r) can be constructed. 
For a deeper understanding of the stationary equation, it is helpful to discuss it 
explicitly. In Appendix B, we study the simplest case, the GUE, as an example. 
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4.4 Diffusion Equation on the Unfolded Energy Scale 

From the previous discussion and from Section 3.4, it is obvious that the transition 
on the unfolded energy is a diffusive process as well. Naturally, the fictitious time 
is now given by 

T = AV2 (4.36) 

and the combination of rescaled energies ^ and source variables j provides the ficti- 
tious coordinate space. We write 

P = ^ + J (4.37) 

with p having the form (3.8) in the case of the GUE and we make the replacement 

^ + J ^ P (4.38) 

where p is of the form (4.13) in the case of the GOE and GSE, respectively. The 
structure of the curved space in which the diffusion takes place remains unchanged 
and so does the Laplace operator, we simply have to replace r with p in Eq. (4.17). 
Hence, the generating functions 

Zf3kip,T) = lim Zf3k(r,t) (4.39) 

obey the diffusion equation 

-AppZpk{p,T) = -T^Zfikip.T) with lim 2;^fc(/>,r) = zf^ip) (4.40) 

for all three symmetry classes and for arbitrary initial conditions. The formal solu- 
tion 

zpkip.r) = exp(^^A^,) z^^^ip) (4.41) 

is completely analogous to Eq. (4.20). In the case of the GUE, a separation of the 
type (4.9) is of course also possible. The diffusive processes on the original and the 
unfolded scale are the same, the initial conditions are different. This is no surprise 
since the Laplacean (4.17) is independent of the level number N which appears only 
in the initial condition. Once more, the space of the eigenvalues of the supermatrices 
proves the appropriate one since it describes just the correlations and does not care 
about the scale on which they are measured. In other words, a once unfolded initial 
condition remains always unfolded, even during and after the propagation into the 
regime dominated by the Gaussian ensembles. It should be mentioned that the 
diffusion equation on the unfolded scale is more general than the one on the original 
scale. This is because A, and thus r, as discussed in Section 3.4, turn out to be the 
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relevant parameter on the unfolded scale for a wide class of parameterizations (2.24) 
on the original scale. 

Since the diffusion equations on the original and on the unfolded scale are the 
same, the integral representation (4.25) has to hold on the unfolded scale as well, 

zpk{p,T) = j Tpk{s,p,T)z^p^{s)Bpk{s)d[s\ (4.42) 

which becomes Eq. (3.32) in the case /3 = 2. 

The stationary equations for the pure ensembles on the unfolded scale can dir- 
ectly be read of from the diffusion (4.40). In the limit r oo, the correlations are 
purely chaotic and thus independent of r. The same must be true for the generating 
function, 

Ymi^zpkip.r) = z^^lip) = ^Imi^ 4^(0 • (4-43) 
Consequently, the derivative with respect to r in Eq. (4.40) vanishes and we conclude 

^Pp4Kp) = 0. (4.44) 

The boundary conditions are such that the correlations are translation invariant, 
this implies that the generating function can depend only on the differences of the 
eigenvalues p. 

By using Dyson's Brownian motion picture [27], hierarchic relations were con- 
structed in Ref. [9] for the correlation functions on the original and the unfolded scale, 
Rk(xi, . . . , Xk,t) and Xk(^i, . . r), respectively, in which the energies provide the 
fictitious coordinates. These equations couple the correlations of the order k to those 
of the order k + 1 and thus cannot be viewed as a diffusion. Due to the application 
of supersymmetry and the graded eigenvalue method, our equations for the gener- 
ating functions describe an independent diffusion for every order k in the extended 
fictitious coordinate space. 

5 Transition from Poisson Regularity to Chaos 

In Section 5.1, we discuss some general aspects of regularity in random matrix 
theory. The Poisson Ensemble is introduced in Section 5.2. The generating function 
is worked out in Section 5.3. In Section 5.4, the correlation functions are evaluated. 
In section 5.5, the two level correlation function is discussed in detail. 

5.1 First Considerations 

The Gaussian Ensembles model quantum systems whose eigenenergies repel each 
other because of some interaction. This famous Wigner repulsion, now often re- 
ferred to as quantum chaos, manifests itself in the simple fact that the matrices H^^^ 
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which form the Gaussian Ensembles are fully occupied. Since no basis is preferred, 
the GOE, GUE and GSE are invariant under 0(iV), C/(iV) and SpiN) rotations, 
respectively. In other words, the probability density functions depend only on the 
eigenvalues. As well known [3], the non-trivial part of the ensemble average becomes 
then an integration over the eigenvalues Xn \ n = 1, . . .,iV involving the Jacobian 
I Ajv(a;^"^^)|'^. The Vandermonde determinant Ajv(a;^"^^) = nra>m(^"-^^ — Xm'^) annihil- 
ates the whole integrand whenever two eigenenergies degenerate. This is the Wigner 
repulsion. The chaotic features of the correlations on the unfolded energy scale are 
almost solely due to this effect. The special form of the probability density function 
is not important as long as it is rotation invariant [3]. This robustness or universality 
of the chaotic fluctuation properties has recently been verified in detailed studies for 
spectral correlations [28] and also for scattering systems [29]. 

How to construct a random matrix model for quantum systems which do not 
have chaotic, but regular fiuctuation properties? This question has been studied 
by numerous authors addressing many different physical situations, especially those 
which show a transition [9] toward chaotic fiuctuation properties. As already stated 
in the introduction, in this work, we are mainly interested in a regular system 
that does not show any spectral correlations. The Poisson Ensemble models such 
a system. In view of the discussion above, it comes as no surprise that ensembles 
without spectral correlations are efficiently designed by breaking the rotation invari- 
ance. Loosely speaking, the infiuence of the Vandermonde determinant becomes the 
weaker the stronger the rotation invariance is broken. This models the vanishing of 
the Wigner repulsion. A transition ensemble of particular interest is the Band Ran- 
dom Matrix Ensemble in which all entries of the matrices beyond a certain distance 
from the diagonal are set equal to zero. For a fixed, large dimension of the matrices, 
the fiuctuations become more and more Poisson regular as the bandwidth is made 
smaller [30, 31]. Closely related and maybe better suited in a wide class of physical 
situations are Gaussian Ensembles which are coupled in such a way that strength is 
transported between next neighbors. Those ensembles consist of tri diagonal matrices 
whose entries are again random matrices drawn from Gaussian Ensembles. Here, 
for a fixed dimension, the ratio between the variances of the Gaussian probability 
density functions for the side and the main diagonal controls the transition. These 
ensembles are being used as models for systems which show localization effects and 
to study the fiuctuation properties of mesoscopic particles [32, 33]. Recently, an al- 
ternative ensemble was studied [34] that is in a certain sense rotation invariant but 
still yields an interpolation between Poisson regular and chaotic fiuctuation prop- 
erties. Here, however, we are interested in transition ensembles of the form (2.14). 
Hence, we have to construct an ensemble of matrices H^'^^ that models pure Poisson 
regular fiuctuation properties. 
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5.2 Poisson Ensemble 



The easiest way to break the rotation invariance is to assign a fixed value to certain 
matrix elements. If we want to eliminate all correlations we simply have to remove 
all interactions from the matrices H^'^^ by setting all offdiagonal matrix elements 
equal to zero, 

N 

Pi°)(i7(°)) = nP^°^(^i°^) n S(ReHi'y)S(lmHi'y) (5.1) 

n=l n>m 

where p(°)(2;) is a smooth, symmetric, but otherwise arbitrary, probability density 
function with the normalization 

+00 

p(°\z)dz = 1 . (5.2) 



This ensemble is referred to as Poisson Ensemble. The lack of any correlation is 
straightforwardly shown. From definition (2.2) we find 



Rf\x) = Np^^\x) (5.3) 



for the level density and 



p=i ^ J* / 

for the A;-level correlation functions. As mentioned already in Section 2.1, we always 
assume Xp Xg for p ^ g in order to remove the trivial contributions containing 
S(xp — Xq). Obviously, there are no true correlations, because the A;-level correlation 
function is simply a A;-fold product of level densities. Notice that the prefactors in 
Eq. (5.4) ensure the normalization of the A;-level correlation function to Nl/(N — k)! 
which holds by definition for arbitrary ensembles. For the same reason, the prefactor 
in front of the product of all level densities is unity in the decomposition of the 
higher order Gaussian correlation functions into the lower order ones according to 
Eq. (2.11). It is therefore not very convenient to define A;-level cluster functions in 
the Poisson case. On the unfolded energy scale, however, all these normalization 
factors become unity and we can write 

Xf)(ei,...,6) = 1 for all A; 

yf)(ei,...,6) = for A; > 1 (5.5) 
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without causing inconsistencies. These equations can also be viewed as the definition 
of a Poisson Ensemble. The inclusion of the real part of the Green function yields 
the contribution of a Cauchy principal value integral, 

+00 

= — [ P (^—] p^°Hz)dz T iNp^^Hx) . (5.6) 

IT J \X — Z J 



-00 



All other formulae above hold the same if we replace Ri(x) with Ri(x). 
5.3 Generating Function 

As the initial condition for the diffusion equation or the the integral representa- 
tion (3.16) we need to evaluate the generating function (3.15) for the Poisson En- 
semble (5.1). We find 



^ ,s) = j P^^\h('^) detg-i (.± ® lA. - hk ® d[H('^] 

Cfl^))"^ (5.7) 

since all diagonal elements are equally distributed. There are various ways to com- 
pute the integral 

-1-00 ^ 

Cr(.) = [p(%)ll'-^^dz, (5.8) 

-00 p=^ ^p^ ^ 

the easiest one is probably to use the identity 

p-l H-p ^ p-i llqj^pKH'q h^'p)h^'P ^ 

which holds for an arbitrary set of variables j-ip^Vp, p = l,...,k. Degeneracies of 
the type /^p = Vq for some p, q are taken care of through cancelations in the right 
hand side. Even degeneracies of the type jip = jiq for some p, q come out correct if 
viewed as the limit jip jiq. The identity (5.9) can be proven by induction or by 
standard methods of complex analysis. Hence, for the integral (5.8) we find 



^ T'^ (0)/ \ ^ 

cf\s) = 1 + Y^bpis) [ p^dz = 1 + ^Y.M')Ri \'p^) (5-10) 

„_i J 'S„-| Z iV 



where the function 

jk 



M^) = -fi - = {^Sp2-Spl) 11 + + (5.11) 

nq^pK^ql ^pl) q^p *gl *pl 
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can be considered as the projection of Bk(s), the square root of the Berezinian, onto 
the coordinate axis Spi. It is interesting to remark that these functions obey the 
sum rule 

k 

J^bpis) = -trgs . (5.12) 
p=i 



Collecting everything we arrive at 



It is easily checked that zj^\x + J) is indeed the generating function for the Poisson 
Ensemble. We have bp(x + J) = 2Jpbp(x + J) with bp(x) = 1, the derivatives of 
bp(x + J) cannot produce expressions proportional to Jp or 1/Jp. Hence, the only 
term that contributes to the correlation functions (2.6) is the one containing the 
product Ylp=i -Jp- Counting the prefactors correctly we find the results given in the 
last subsection. 

Before we can compute the generating function (3.33) on the unfolded scale, 
we have to specify the one point function R^^\x)^ in particular its level number 
dependence. As discussed in Section 2.2 the two level densities involved, R^^\x) 
and R^^\x)^ ought to have the same asymptotic N behavior. This suggests to 
choose the GUE level density for the Poisson Ensemble. We set R^^\x) = R^j^\x) 
which implies through the Kramers Kronig relation also the equality of the real 
parts. For large N , we can use the Pastur equation of Section 3.3 to determine the 
one point function for the transition which is, as discussed, just given by Eq. (3.29). 
On the unfolded scale, we are only interested in the transition parameter a of the 
order l/-\/iV". Consequently, the a dependence vanishes for large N and we have 

Ri{x,a) — > Rf \x) = R^^\x) — > - (x =F iV^N - xA (5.14) 

and thus D = 1/Ri(0,a) = 7r/-\/2iV". The equations (5.7) and (5.10) yield in for- 
mula (3.33), 

4°)(5) = exp Qirn^ iV In Cr(i?5)) = exp ||^lim^ iV In ||l + ^^ ^ T&,(5) j j 

k 

= Y[exp(^iTTbp(s)) (5.15) 

p=i 

where we used bp(Ds) = Dbp(s) and R^f^Dspi) ^i/D for N oo. The signs 
are determined by the choice of the sign of the imaginary increment in the Green 
function. 
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5.4 Correlation Functions 



In order to evaluate the A;-level correlation functions, we insert the result (5.13) into 
Eq. (3.18) and arrive at 

(-1)*^ [ 



k \^ 



(^1 + ]^E^?'(W(v)j B,{s)d[s] (5.16) 

for non-zero a. As pointed out in Section 3.2, the limit a ^ yields the correct 
result for the generating function, i.e. here in this case the generating function (5.13). 
As a support of the general discussion in previous sections and as a further test 
of our calculation, it is worthwhile to check the opposite limit for the correlation 
functions (5.16). Following Section 2.2 we are free to replace the eigenvalues s in the 
initial condition with s/Vl — a^- Since Z^\s) is a continuous function oiQ^\s)^ we 
only need to check the limit in the rescaled function C^\sl\/l — a^). It is instructive 
to do that explicitly in the form (5.10). Because of the homogeneity of h,p{s) the 
relevant limit reduces to 



lim Ri , ' = lim — / , , ^dz = r (5.17) 

and gives after resummation of the expansion in Eq. (5.10) the expression 

k 

AO) ' 




lim Zr ( ^== ) = n ( ^ 1 (5.18) 



which is indeed the correct result [18] for the GUE. Notice that the limit (5.17) 
is just the unaveraged Green function (3.26) which has already been used in the 
discussion of the Pastur equation in Section 3.3. 

Lenz [11] presented a calculation of the correlation functions for the Poisson GUE 
transition employing the Mehta Dyson method, i.e. without using supersymmetry. 
After several changes of variables, he found an integral representation for the one 
level correlations which coincides exactly with our formula (5.16) for k = 1. For 
the two level correlations, he constructed an integral representation which contains 
a ratio of level number dependent determinants. Unfortunately, this representa- 
tion seems not to be amenable to further analytical treatment. In particular, it is 
not clear how to take the limit of infinetly many levels. Hence, the advantage of 
supersymmetry and the graded eigenvalue becomes obvious in the simple form of 
the generating function (5.13). The crucial quantities are the functions bp(s) which 
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are related to the Berezinian. They generate the determinants in Lenz's integral 
representation and allow to take the limit N ^ oo directly in the form (5.15). 

The construction of the physically interesting correlations Rk(xi, . . . ,Xk,a)is not 
completely straightforward. In order to apply the procedure described in Section 2.1, 
the contribution from the Poisson Ensemble in the integrand has to be expanded 
in a binomial series. No further simplifications seem possible in this procedure. 
Fortunately, things become simpler on the more important unfolded energy scale. 
For non-zero A we find from Eqs. (3.35) and (5.15) 

Xfc(ei,...,6,A) = J Gk{s-C,\)llexp{^nrbp{s)) Bk{s)d[s] (5.19) 

p=i 

which implies certain simplifications due to the product structure in the integrand. 
Using the symbol we can write 

XkiCi, • • • , 6, A) = J Gk{s - e, A) S n exp{^nrbp{s)) Bk{s)d[s] (5.20) 

p=i 

representing the A;-level correlation functions on the unfolded energy scale for the 
transition from Poisson regularity to chaos as a 2A;-fold integral. Since the generating 
function (5.15) depends only on the differences of the eigenvalues, the correlation 
functions are indeed translation invariant as discussed in Section 3.4. 

In the case of the level density, it is very easy to see that the initial condition 
takes the form 

?y2;|°^(s) = sin7r5i(s) = sin 7r(isi2 — sn) (5.21) 

and is thus according to Appendix C of the form that ensures the validity of 
Eq. (3.40) such that the level density is unity everywhere for all values of A as 
it should be. 

For k > 1 the functions bp(s) have differences of eigenvalues Spi in the denom- 
inator which still carry imaginary increments. The exponentials of these functions 
have to be viewed as operators that generate distributions which complicates the 
integrals to be performed. However, the two level correlations can still be expressed 
as a double integral. 

5.5 Two Level Correlations 

For the analysis of data, the two level correlation function is of particular import- 
ance. The general considerations of Section 3.5 apply and the function X2(r, A) 
with r = ,^1 — ,^2 is of the form (3.42). Thus, we have to evaluate ^Z2^\o,ti,t2), 
i.e. S exp(=pi7r5i(s)) exp(=pi7r52(s)) in the coordinates (3.41). In Appendix E we show 
that the initial condition takes the form 

S4°)(0,ii,i2) = ^Ke (^exp - 1^ (5.22) 
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where the exponential has to be interpreted as a power series involving the operator 
Using the principal value and the ^-function, the initial condition can be 
written as 

^z, (U,ti,t2j - 2 A. 22-(2n)!(2n-l)! dtl^-' \h 



+ 2^^ 22-+i(2n+l)!(2n)! ^ ^^^ ^^ 

which is obviously not very convenient. We therefore express the initial condi- 
tion (5.22) by means of an integral transform 



°" ^-l)-(7r(i2 +i2))2n+l Q2n 



(o)/n ^ ^\ _ ^2) ( /o_/^2 , ^2\„A exp(-eM) cos(iiM) 



34^^(0,^1,^2) = ^ 4 h{j2Tr{tl + tl)u\ ^ '^^ ' du 



(5.24) 

where Ii(z) is the modified Bessel function of first order. The convergence is ensured 
by the infinitesimal increment e. This representation allows, inserted into Eq. (3.42), 
the evaluation of one of the three integrals, details are given in Appendix F. We arrive 
at the double integral 

00 00 / ^ X 

4 ^ /• , / 1 1\ /■ , I P \ 



^^^^'^^ = ^^°^y^"U + ^2r/A2+^jy^^^"PV 2A2/ 
^ ^ 

/i (\/'2Trup) / , \ 

— ^ J2 (^V^2 + i2Mr/A2^j (5.25) 

in which J2iz) is the Bessel function of second order. Since the u integral converges 
as it stands if the p integration is done first, we have already taken the limit e ^ 0. 
Thus, for the first time, a closed formula for the two point correlation is given which 
is valid for all values of the transition parameter A. Yet another double integral 
representation can be found 

00 

X2(r, A) = — Im J du exp (it \^u) ~ 


I dr]!]"^ exp I- (^u - i2-^^ 

^ ^ 

Ji (^V2^\^u^/u - i2r/A2?7) 



/2^A2mVm - i2r/A2 



(5.26) 
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which involves only one Bessel function. The derivation is sketched in Appendix F. 
According to the integral representations, the two level correlation function has the 
property 



X2(0,A) 



= (5.27) 



which simply means that an arbitrarily small GUE admixture to the Poisson spec- 
trum lifts immediately all degeneracies. 

Our result for the two level correlation function can be checked by discussing 
some limit relations. First, for large energy differences, the two level correlations 
must approach unity, 

lim X2(r,\) = 1 . (5.28) 

r^oo 

Second, the limit of vanishing transition parameter must also yield unity 

limX2(r,A) = 1 (5.29) 

since the Poisson spectrum is totally uncorrelated. It is shown in Appendix G that 
our result indeed satisfies both limit relations. Furthermore, it is important to 
compare our formula for small but finite values of the transition parameter to the 
perturbative calculations of Refs. [9, 10]. From Eqs. (3.45), (3.46) and (5.24) we find 
for small A 

X2(r,A) ~ X2,o(r/A) = L J ^xp 1-^] sin y rffc (5.30) 

^ ^ 

in perfect agreement with the results of Refs. [9, 10]. Notice that the definitions of 
the transition parameter differ, A in Ref. [10] has to be replaced by A^/2 to compare 
the formulae. Details of the derivation of Eq. (5.30) are given in Appendix G. 

Since further analytical treatment of the integrals (5.25) and (5.26) seems not 
possible, we resort to a numerical calculation of formula (5.25) for some values of 
the transition parameter A. In Fig. 1 we show the results for the three A values 
0.1, 0.5 and 0.7. The numerics is not completely trivial, the p integration yields 
an oscillating function with large amplitudes as integrand for the u integration. 
Therefore, for the larger A values 0.5 and 0.7, the calculation is only meaningful up 
to energy differences of r = 1.2. 

For A = 0.1, the uncorrelated Poisson spectrum has been affected by a rather 
weak level repulsion. All degeneracies are lifted, but the correlations are only relev- 
ant for smaller energy differences r. The function has therefore a very steep gradient 
near the origin r = 0. However, the level repulsion has to be compensated for in- 
termediate r before the Poisson situation is restored. This compensation manifests 
itself in in a region where A2(r, A) exceeds unity which we will refer to as the over- 
shoot. This effect has already been observed in Refs. [9, 10]. It can be understood 
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Figure 1: Two level correlation function of the transition from Poisson to GUE 
for different values of the transition parameter A, calculated numerically using for- 
mula (5.25). In the Poisson case, i.e. for A = 0, the function is unity which is 
not drawn. The three thin solid lines correspond from the left to the right to A 
values of 0.1, 0.5 and 0.7, respectively. The thick solid line is the pure GUE case 
corresponding to A ^ oo. 
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as a consequence of a normalization or sum rule. In the Poisson case, the integral 
over the two level cluster function on the unfolded is zero which implies that the 
integral of Y2(r, A) = 1 — X2(r, A) is approximately zero for small A thus causing the 
overshoot. In fact, it is easily shown that the approximation (5.30) fulfils the sum 
rule, i.e. the integral of 1 — X2fi(r/X) yields zero. On the other hand, the integral 
over the two level cluster function on the unfolded scale is unity in the GUE case. 
Hence, the integral of Y2(r,X) changes from zero to one as the transition is made 
from A = to A ^ 00. We mention in passing that on the original energy scale 
where the level number N is finite, the two level cluster function must satisfy the 
same sum rule for all values of the transition parameter due to the definition (2.16). 
Thus, the sum rule and the limit of infinitely many levels iV ^ oo do not commute 
in the Poisson case. However, the interpretation of our results in terms of a change 
of the sum rule on the unfolded energy scale is very useful. It implies immediately 
that the overshoot for smaller A has to disappear as A is made bigger. Indeed, Fig. 1 
shows that X2(r, A) for A = 0.5 exceeds unity only very slightly. For A = 0.7, no 
overshoot is visible anymore. 

Furthermore, it becomes obvious from Fig. 1 that the two level correlations reach 
the GUE limit rather quickly when the transition parameter A approaches values of 
about unity. This is consistent with the numerical simulations of this transition in 
Ref. [7]. Roughly the same speed has also been observed and analytically shown for 
other, but related transitions, in particular for the breaking of time reversal invari- 
ance, i.e. the GOE GUE transition [12], and the breaking of quantum numbers [15]. 

We mention finally that double integral representations of the level number 
variance S^(i,A) for the transition from Poisson regularity to GUE can be found 
straightforwardly by using the general expression (3.48) and the methods of Ap- 
pendix F. Since the results are similar to formulae (5.25) and (5.26) but more com- 
plicated, we decided not to present them here. 

6 Summary and Discussion 

We analyzed the statistical model of a system that can be written as the sum of 
an arbitrary and a chaotic part. This choice is motivated by a whole variety of 
experimental and theoretical investigations. The chaotic part was drawn from the 
Gaussian Ensembles whose fiuctuation properties are known to describe chaos in 
quantum mechanics. We expressed the generating function of the spectral correla- 
tions for the transition from the arbitrary to the chaotic regime as a supersymmetric 
matrix-model. A Pastur equation for the level density involving the level density 
of the arbitrary part was found through a saddlepoint approximation. In the case 
of broken time reversal invariance, the graded eigenvalue method was used to de- 
rive a closed integral representation of the generating function which contains the 
generating function of the arbitrary correlations in the integrand. The number of 
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integrals to be done is twice the order of the correlation. The unfolding procedure 
could be performed directly in this integral representation which immediately im- 
plied the translation invariance of all correlations, provided the correlations of the 
arbitrary part were translation invariant. Employing this insight, we showed that 
the level density is, as expected, unity for all values of the transition parameter. 
Furthermore, the two level correlation function could be written as a double integ- 
ral. Consequently, a double integral representation for the level number variance 
was found, too. 

For all three universality classes, we derived diffusion equations for the generat- 
ing functions which are completely equivalent to the integral representations. The 
diffusion takes place in the curved space of the eigenvalues of the supermatrices. 
Only the radial part of the Laplacean contributes such that no anticommuting vari- 
ables have to be considered explicitly. The squared transition parameter provides 
the diffusion time. The generating function of the arbitrary correlations is the initial 
condition. Hence, arbitrary correlations are propagated diffusively into the chaotic 
regime. The form of the diffusion equation is identical on the original and the un- 
folded scale which is a strong demonstration for the utility of the graded eigenvalue 
method. We expect that these diffusion equations are also valid if the initial con- 
dition is not purely spectral, for example, in the case of a coupling to channels. 
Stationary equations for the pure chaotic correlations follow directly from the dif- 
fusive process. 

We discussed the transition from Poisson regularity to chaos by taking an initial 
condition without spectral correlations. The generating function can be worked 
out in closed form and was used in the general integral representation. On the 
unfolded energy scale, the generating function is an infinite series of distributions 
which could be written as an integral transform. Using the general formula, the two 
level correlations were expressed as a double integral. This is the first time that such 
an integral representation has been found. Most likely. There is no other method 
that allows for an evaluation of both ensemble averages and yields a low dimensional 
integral representation at the same time for all values of the transition parameter. 

Appendix 

In Appendix A, the Pastur equation is derived. The stationary equation for the 
GUE is discussed in Appendix B. In Appendix C, the level density on the unfolded 
scale is evaluated. The two level correlations on the unfolded scale are discussed in 
Appendix D. In Appendix E, the Poisson initial condition for the two level correl- 
ations is worked out. Integral representations of the two level correlation function 
are constructed in Appendix F. In Appendix G, checks of the two level correlation 
function are performed. 
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A Derivation of the Pastur Equation 

Since the generating function z[^\s) is independent of the superunitary matrix u, 
we need only the radial part of the gradient operator d/da. In Ref. [35] the full 
operator was expressed in spherical coordinates, the radial part reads 



1 / 9 . d \ 

U \ Til- h iT22Tr- U 

OSi dS2j 



(A.l) 



where the 2x2 matrices T,p,p have unity in the entry pp and zeros elsewhere. We 
omit the index k = 1. From Eq. (3.21) we find that the saddlepoint (Tq = Uq^sqUq 
is given by 







trg ^(so 



uodauQ ^ . (A. 2) 



Since the supertrace has to vanish for all matrix differentials da, the diagonal matrix 
in front of uodau^^ is zero and we obtain the two coupled equations 











2 d 
■ — (sio-a;) + — In Z[°\s) 
a'^ osi 



S = So 



2 d 
— (is2o - x) + i— In Z\ 



a 



(A.3) 



S = So 



By setting si = x' — J' and is2 = x' + J' we introduce new variables x' and /' 
which play formally the role of energy and source variable in the generating function 



Z'l 



(0), 



The difference of the Eqs. (A.3) yields 







-^{sio - is2o) + -^InZl 



a 



(A.4) 



S = So 



which is solved by sio = is2o, i.e. /' = 0, because the derivative of the generating 
function with respect to the energy is zero for vanishing source term. The sum of 
the Eqs. (A.3), 



= H — ^(sio + is20 



(A.5) 



S = So 



involves now the derivative of the generating function with respect to the source 
variable at /' = which is according to Eq. (2.6) nothing else but the level density. 
This gives Eq. (3.22). 

In order to integrate over the massive modes we have to expand up to second 
order, 

1 d'^L 

L(ao + da, x + J,a) = L(ao, x + J,a) + — trg- 



'da^ 



da , 



(A.6) 



(7 = (7o 
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and evaluate the Gaussian integral over da. It is straightforward to work out the 
second order term explicitly. However, this is not necessary since the gradient is 
Hermitean and the Lagrange function is well behaved. Thus, the matrix /da^ at 
(7 = (To is also Hermitean and can be absorbed into da. Consequently, the Gaussian 
integral converges and yields just unity. The generating function 

Zi(x + J,a) = exp L(ao, X + J, a) 

= exp (^-^trgiso- X - Jf^ (N ^ oo) (A.7) 

gives the result (3.23). 



B Stationary Equation for the GUE 

The ansatz (4.9) reduces the diffusion (4.8) in a curved space to a diffusion in a 



Cartesean space. Similarly, the ansatz Z^j^^r) = W2]} (r)/ Bk(r) reduces Eq. (4.33) 



to the Cartesean stationary equation 



+ W^lHr) = (B.l) 



where we used that r -^In Bk(r)/dr is equal to the constant —k. An equivalent form 
of this equation is given by 



Q^ + r) W2k(r) = r'W^,'(r) (B.2) 

whose left hand side suggests the ansatz W2]}{i^) = 6xp( — r^/2)F2^^^(^)- It leads to 
the equation 

- r') yil\r) = (B.3) 

which is of oscillator type. Obviously, the eigenvalues r^i and irq2 can be arbitrarily 
paired such that the solution can be built with functions f ^^^(r^i, irg2) that satisfy 

^ rli] v^'Hrpi.^r,2) = f 1^7^ - (^^2)^') v('Hr,u^r,2) . (B.4) 



\(J^v\ ) \d(trq2) 

The equation is separable and the solution is a product. The addition of the term 
(2ra+l)f (^)(rpi, irq2) to the equation makes both sides to the oscillator wave equation 
of index n. The imaginary unit in front of irq2 takes care of a proper cancelation of 
the ra-dependent parts. Thus we find for arbitrary n the solution 

v^^Hrpi,irq2) = ipn(rpi) ^n(irq2) (B.5) 
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where the functions (p^ are a linear combination of both fundamental solutions of 
the oscillator equation [25]. Every linear combination of these product solutions is 
again a solution, in particular we can choose 

JV-l 

■v^^Hrpi,irq2) = KN(rpi,irg2) = ^ ^nl^-pi) ¥'n(«?'g2) (B.6) 

n=0 

for arbitrary N . Moreover, all boundary conditions are satisfied by 

^ifcV) = det[li'jv(rpi,irg2)]p,g=i,...,fc (B.7) 

which is indeed the correct result. Since none of the equations involved is expli- 
citly iV-dependent, one can also construct envelope solutions corresponding to an 
asymptotic expansion in N . 

C Level Density on the Unfolded Scale 

By construction, the level density of the arbitrary ensemble has to be unity on the 
unfolded scale, x['^\^) = 1. According to Eq. (2.6), this is only possible if the 
generating function has a power series expansion of the form 

oo 

CJ^(0)(5) = - tS2) + J2 ''nisi - tS2r (CI) 

n=2 

with some coefficients c„. We drop the index A; = 1 in our notation. For the level 
density of the transition, we find from Eq. (3.35) the double integral 

+00+00 



Si - IS2 



00 — 00 
+ 00 + 00 

2 / i-Ti [st + s^2 



1 ; ; / ^'2,2 



ttX^ J J ' \ \ 

— 00 — 00 



-j^ 00 



l--J2''n(si-tS2r-^] ds^dS2 . (C.2) 



n=2 



Apart from the Gaussian functions, the integrand depends only on the difference 
si — is2, the sum si + is2 does not appear. Thus, we can apply the formula 



+00+00 

1 

exp 



7rA2 

— 00 — 00 



(^-^{sj + sl)^ f{Si-tS2)dSidS2 = /(O) (C.3) 



which holds for any well behaved function /. This is a well known mean value 
theorem from complex analysis. It can be shown by introducing polar coordinates 
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Si + is2 = uexpii-d) and then expanding the integrand in the phase. In the limit 
A ^ 0, it is also valid due to the occurrence of S(si)S(s2). Formula (C.3) yields 
immediately the result (3.40). 

D Two Level Correlations on the Unfolded Scale 

With the new integration variables (3.41), the relevant parts of the integrand read 

trg(s-0^ = ^{sj + si + tj + tl) - r{si - is2) - r{h - it2) 



(si - iS2 +ti - it2){si - iS2 + ^1 + it2) 

4it2 



(si - is2 - h - it2){si - is2 - h + it2) 



.(D.l) 



As in the case of the level density in Appendix C, apart from the Gaussian functions, 
the whole integrand depends only on the difference si — is2. We can thus use 
formula (C.3) and find 



+ 00 + 00 



A2(r,A) = l I exp^-^itl + 4)-^ih-tt2: 



ii it 



1"2 (0)/ 



r.2^,2^2 ^^r(0,^i,^2)rfiirfi2 (D.2) 
I'l + '■2) 

which does not contain r anymore as required. Since the integrand is odd in each of 
the variables ti and ^2 except of the r dependent exponentials, it is enough to keep 
the odd parts of those which gives formula (3.42). 

E Poisson Initial Condition for the Two Level Correlations 

In the case k = 2 we have to consider the functions 

, .X {iSl2 - Sll){iS22 - Sii) {iSl2 - S2l){iS22 - S21) 

bi{s) = _ ± and 62(5) = ± _ ± 

•^21 •^11 •^11 •^21 

(E.l) 

in the exponent. We introduce the short hand notation indicating the sign of 
the imaginary increments. From Eq. (2.5) we find the expression 

(i2)2s4°^(s) = exp(+i7r(5^- +5^-)) + exp(-i7r(5++ + 5++)) 

-exp(-i7r(5+- - 5+")) - exp(+i7r(5-+ - 5"+)) . (E.2) 
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which can be simplified with the sum rule (5.12) giving 



(s) = —Re ^exp(i7r(5-^ — ^2 "*")) ~ exp( — ivrtrgs)^ 



2 V / 2 



= — Re exp( — ivrtrgs) ^exp( — i27r52 "*") — 1) • (E-3) 

We use the coordinates (3.41) and put si — is2 = 0. The supertrace vanishes due 
to trgs = si — is2 and we arrive at Eq. (5.22). Obviously, we have to deal with a 
function of the operator which is defined through its power series expansion. 

By expanding the exponential and using 

—— = ) ' -— = \ ' .^^ ^ p ± +ivr^(ii) (E.4) 

we arrive at Eq. (5.23). 

F Integral Representations of the Two Level Correlations 

The integral (5.24) is a special case of formula 6.643.2 given in Ref. [37]. It can also 
be derived directly by expanding the exponential in Eq. (5.22), using Eq. (E.4) and 
the formula 

oo 

— = i J exp( — it^u)du . (F.l) 
^ 

Equation (5.24) is inserted into the integral representation (3.42) and polar coordin- 
ates ti = pcos-d, t2 = ps'm-d are introduced. The -d integration is of the form 

exp(ia cos + i5 sin f^) cos sin rff^ = — 27r ^ ^, Jq i\/ + 6^ ) 

narlh V / 



dadh 





27r 



-^^.J2[Va^) (F.2) 



where Jn(z) is the Bessel function of order n. In our case we have a = —p(u — ir/X^ ) 
and b = pr/X^. Collecting everything yields the integral representation (5.25). The 
p integration which has the form 

oo 

M(a,5,c) = j exp(-cp'^)Ii(ap)J2(bp)dp (F.3) 



can be rewritten by making use of Sonine's formula [38] 

J^+.+i{z) = ^^^^^^ ^ i J^{z sm{}) sm''+'{}cos^'+'{}d{} (F.4) 



44 



for 11=1 and z/ = 0, 

1 

,2 



.J^(z) = z j .h(zri) ri' dri (F.5) 



where we changed the integration variable to 77 = sin The order of the Bessel 

function J2{bp) is thus reduced by one which allows the application of Weber's 
formula [38] 

I exp(-cp^)I^(ap)J^(bp)pdp = ^ exp | " ^ j J<y (^^^ (F.6) 

^ ^ 

for z/ = 1. Collecting everything we find 

M(„.6.c) = Ae.p (h!) /exp (-^) .7, „) ,F.7) 



which yields with a = V2vrM, 5 = — i2ur/\'^ and c = 1/2A^ the integral repres- 
entation (5.26). 

G Checks of the Two Level Correlations 

We start with a general expression for the two level correlation function, but instead 
of Eq. (3.42) we use the equivalent form (D.2) and shift the integration variables 
according to ^ + r and ^2 ^ ^2 — ii" which gives 

+00+00 

X. / ^^ f f 1 /.2 .2^^ (h - r'Kh + ir) 

-00 — 00 

'^z^^\0,ti - r,t2 + ir)dtidt2 . (G.l) 



The contribution from the Berezinian yields i/4(ii — 7^2)^ in the limit of large energy 
differences r. For a Poisson initial condition, we have from Eq. (5.22) 

s4°V,ii-r,i2 + ir) = iRe|^exp|^-i|(ii-ii2)(^l + i^^^^^ -1^ (G.2) 

where now, for consistency reasons, the imaginary unit in front of ^2 which appeared 
due to the shift has to be ignored when taking the real part. Hence, the limit of 
large r yields 

1 TT 

-Re (exp (-7(^1 - 7^2)) - 1) = - sin^ -(^1 - 7^2) • (G.3) 
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Collecting everything we find 

+00+00 



lim X2(r,A) 



1 



7r3A2 



exp 



2A 



(G.4) 

which gives the limit relation (5.28) by applying the integral formula (C.3). 

The representation (G.l) of the two level correlation function allows to consider 
the limit of vanishing transition parameter, too. For A ^ the Gaussians become 
the functions ^(^1)^(^2)- After integration we find 



limX2(.,A) = 4 lim - ^^1^ 



1 



(G.5) 



and thus the limit relation (5.29). 

Finally, we discuss the two level correlation function for small but finite values 
of the transition parameter. From Eq. (5.24) we have to zeroth order 



A2=0 



16 



(G.6) 



which has to be inserted into formula (3.46) for ra = 0. Luckily, there are some 
cancelation that decouple the two integrations to be performed. We arrive at 
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M^iri) 
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t2 sin 77^2 dt2 . 
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Integration by parts yields three ordinary principal value integrals for Mi (77), two 
of them cancel each other and we find 



Ml (77) 



1 



+ 00 



/27r J """^^y 2 j h 

—00 ^ ' 

while the second integral gives trivially 

M2(i]) = 77 exp 



^1 \ sinh 77^1 



dti 



?7_ 

' 2 



(G.9) 
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Thus, the product contains Dawson's integral [37], 



^2,o(r/A) = 




(G.IO) 



which is identical to the result (5.30). 
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